High-pressure molecular phases of solid carbon dioxide 
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We present a theoretical study of solid carbon dioxide (CO2) up to 50 GPa and 1500 K using first- 
principles calculations. In this pressure-temperature range, interpretations of recent experiments 
have suggested the existence of CO2 phases which are intermediate between molecular and covalent- 
bonded solids. We reexamine the concept of intermediate phases in the CO2 phase diagram and 
propose instead molecular structures, which provide an excellent agreement with measurements. 

PACS numbers: 62.50.+p, 64.70.Kb, 61.50.Ah, 78.30.-j 



Understanding the evolution of the bonding properties 
of molecular crystals as a function of pressure is a funda- 
mental question in condensed matter physics. Recently, 
this question has received a widespread attention due to 
key progress in experimental techniques and hence, the 
availability of new data for a number of molecular solids 
lj. In the case of CO2, the interest has been further in- 
tensified due to its importance for planetary science and 
technical applications. However, despite several exper- 
imental investigations of the CO2 high-pressure phases, 
the changes of intra- and inter-molecular bonds, as well 
as of electronic and vibrational properties, as a function 
of pressure are not yet well understood. In addition, large 
portions of the CO2 phase diagram are unexplored from 
a theoretical standpoint. 

At low pressure (P) and temperature (T), CO2 con- 
denses as a molecular solid in the cubic Pa3 structure 
(known as dry ice, phase I), which is characterized by 
strong double bonds (C=0 distance of 1.16 A) and rather 
weak inter-molecular interactions . Recently, a high-P 
and high-Tphase of CO2 (phase V) has been discov- 
ered U 0, IH ' which is completely different from 
the molecular solid. It has a polymeric quartz-like struc- 
ture, and a very low compressibility (the experimentally 
derived bulk modulus, B , is 365 GPa). Following this 
exciting discovery, two additional new phases have been 
reported d, El E3 experimentally (II and IV in Fig.^l, 
both of which were credited with unusual properties; 
their structures were described as intermediate between 
that of molecular and covalently bonded crystals. 

In this Letter, we present a study of the CO2 phase 
diagram based on first-principles density functional the- 
ory (DFT) calculations. Our results challenge the inter- 
pretation of CO2-II, III, and IV as exhibiting dramatic 
differences in the nature of the molecular bonding with 
respect to the low-pressure molecular crystal. In place 
of the previously proposed structures for CO2-II and IV, 
we suggest new, molecular ones. We demonstrate that in 
addition to being stable, our newly proposed structures 
give results in excellent agreement with measurements 
and provide a consistent explanation of experimental ob- 
servations. 

At ambient T, CO2-I undergoes a pressure-induced 



transformation to the orthorhombic Cmca symm etry 
(phase III) between 12 and 22 GPa |1 El ES El El El 
Il6| . The exact nature of this transition and the structure 
of phase III are still unsettled; experimental studies have 
suggested a region of co-existence between CO2-I and 
III [lj] , and also an intermediate distorted low-T phase 
flfij |. Furthermore, it was reported that the I to III 
transition is strongly temperature dependent, becoming 
abrupt above 400 K, and that CO2-III may actually be 
metastable. The two newly discovered phases II and IV 
are believed to be thermodynamically stable and quench- 
able. It was proposed that CO2-II has the P^/mnm 
symmetry with elongated C=0 bonds, 1.33 A, reduced 
inter- molecular distances below 2.38 A (hence named as- 
sociated or dimeric phase), and Bq = 131 GPa jfll. floj|. 
CO2-IV was described as Pbcn, with C=0 bonds as long 
as 1.53 A, and bent (0=C=0 angle of 160°) and strongly 
(dipole) interacting molecules || . Such remarkable prop- 
erties would imply that CO2 is losing its molecular char- 
acter in a "gradual" way, a concept promoted as a key for 
understanding its phase diagram in analogy with solids 
like N2O and Si02- In addition, CO2-III, despite having 
a stable C=0 bond, was also assigned a Bq as high as 
87 GPa 9] . The stability of the proposed bent and asso- 
ciated phases, and the unusual bulk properties of CO2-II, 
III, and IV have not been investigated theoretically up 
to date. 
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FIG. 1: Phase diagram of CO2 according to Ref 0, where 
CO2-III (Cmca) is reported to be metastable. 
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We performed a series of first-principles calculations, 
including full structural optimizations, phonon spectra, 
and free energies, in order to study the stability and 
properties of the phases proposed experimentally up to 
50 GPa and 1500 K. The DFT calculations were car- 
ried out within the Perdew-Burke-Ernzerhof (PBE) 01 
generalized gradient approximation (GGA) using the 
ABINIT code which implements plane- wave basis 
sets [lflj . 

First, we have examined the stability of the bent and 
dimeric structures for phases II and IV. Starting from the 
previously proposed PA 2 /mnm and Pbcn with elongated 
molecular bonds, we carried full structural optimizations 
at various pressures up to 50 GPa. In all cases, upon 
relaxing the atomic coordinates the C=0 bond lengths 
decreased by about 15% and 30% for PA 2 /mnm and 
Pbcn respectively, to become comparable to the free CO2 
molecule [2jj. In the case of Pbcn, the 0=C=0 angle 
also straightened to 180°. The energy differences between 
the theoretically stable molecular structures and the as- 
sociated and bent phases are respectively more than 3 
and 6 eV per molecule. These are energy scales corre- 
sponding to the breaking of a covalent bond. They are 
beyond the errors of the GGA and present a strong ev- 
idence that the previously proposed bent and dimeric 
structures are not stable. In order to further investigate 
this issue, we examine below published experimental data 
for CO2-II, III, and IV, and demonstrate that they can 
be explained in terms of two stable structures, Cmca and 
PA 2 /mnm, which are strictly molecular. 

Computed equation of state (EOS) at T — for the 
three molecular structures Pa3, Cmca, and PA?/ mnm 
are reported in Fig. along with the available exper- 
imental data. In all three cases, the calculated EOS 
agrees very well with the measured one; the computed 
PV curves fall within the experimental uncertainties. We 
find about 2% relative volume reduction from Pa3 to 
Cmca 21] above 30 GPa and another 0.5% from Cmca 
to P4 2 /mnm (see insert in Fig. [2}. Yoo et al. reported 
a ~ 5-7% volume decrease associated with the III to II 
transition; however, they also observed lattice strain in 
their phase III samples 0, which may explain the dif- 
ference between their results and ours. Furthermore, the 
PV curve of Cmca comes close to and eventually merges 
with that of Pa3 below 20 GPa in agreement with Aoki 
et al. [14). Thus, the measured EOS data are reproduced 
well despite the lack of molecular association in our op- 
timized PA 2 /mnm. The computed EOS parameters are 
summarized in Table [i] together with values derived form 
experiment The agreement is good for Pa3, but 

there is an order of magnitude difference for the values 
of Bo of Cmca and PA? I mnm. Since we are using the 
same EOS fit as Ref. lid and we in fact agree very well 
with the direct PV measurements, the large difference 
likely comes from extrapolating the experimental data, 
spread above 20 GPa, down to ambient pressure. 
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FIG. 2: Pressure-volume dependence of selected CO2 struc- 
tures: Pa3 (black), Cmca (blue), and P^/mnm (red). The 
solid lines are ab initio calculations; the points indicate ex- 
perimental data from Ref. [Til (no error bars), Ref.0 (Pa3 and 
Cmca), and Ref. HH (Pfa/mnm). The inset shows calculated 
reduction in volume of PAz/mnm and Cmca relative to Pa3. 

The published experimental evidence invoked to pro- 
pose that CO2-II exhibits the properties of a dimeric 
polymorph consists of: (i) a large splitting of the internal 
symmetric stretching mode, v\ [9j, (ii) a broad librational 
mode identified as Bi s |9LJrOJ, and (Hi) powder x-ray 
diffraction measurements llOj. We carried out calcula- 
tions of the vibron spectrum [23| of molecular PA 2 /mnm 
and Cmca as a function of P. A plot of v\ in Fig. 
shows that the measured splitting, as well as the rel- 
ative values for PA 2 /mnm and Cmca, are reproduced 
remarkably well by the theoretical structures. Yoo et 
al. assumed that the large vibron splitting observed in 
CO2-II is an evidence for decreased inter-molecular dis- 
tances; our results indicate that the crystal field in the 
molecular PA 2 /mnm is sufficient to explain the splitting. 
One should note that phonons are calculated as a second 
derivative of the energy; therefore the agreement with 
experimental frequencies within meV, as found here, is a 
strong indication that computed total energies and forces 
are extremely accurate. 

The computed pressure dependence of the Raman- 



TABLE I: Equation of state parameters for selected CO2 
phases. The theoretical data is well fitted to a third-order 
Birch-Murnaghan equation of state. 
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FIG. 3: Computed Raman-active modes (solid symbols) of 
the theoretically stable structures, (a) Symmetric stretching 
vibron. The dashed lines indicate experimental data from 
Ref. [j| they are from a Fermi resonance band, and have been 
shifted by 52 cm -1 for comparison, (b) External modes. 



active external modes of the three theoretical structures 
is shown in Fig. [3t>- The agreement with experiment fl6| 
is again good, though our values are consistently slightly 
lower than the measured ones. The experimental fre- 
quencies of the two Raman-active modes of CO2-II at 
19 GPa are about 260 and 320 cm -1 - the former being 
a broad pea k - and were classified as B l£ , and E g respec- 
tively |9L llO^I . The broad mode was previously associated 
with dynamical disorder in the lattice. The molecular 
P^/mnm has modes with frequencies 245 and 300 cm -1 
at this P; however it is the degenerate E s which has the 
lower frequency, and therefore corresponds to the experi- 
mentally observed broad peak. Its broadening can be ex- 
plained by invoking a lifting of the E s degeneracy when 
the tetragonal cell is deformed into an orthorhombic one 
(resulting in the Pnnm symmetry). Such a distortion is 
consistent with the x-ray data reported in Ref. |lfj and we 
estimate that between 20 and 30 GPa there is 10 cm -1 
splitting for every 1% modifiaction of a and b. 

We now turn our attention to the x-ray analysis. In 
P4: 2 /mnm, the atomic positions are C(2a) at [0,0,0], 
and 0(4f) at [a;,a;,0]. The associated phase, according 
to Ref. 10, corresponds to x = 0.2732, while our stable 
structure corresponds to x — 0.23075, i.e. x w 0.25 ± 8 
for both structures, where 8 = 0.02. Among the observed 
diffraction peaks, only the (101) and (211) reflections de- 
pend on the sign of 8; their measured intensities relative 
to the calculated intensities of the same reflections are in 
ratios of 2.3 and 0.3 for the associated, and 0.4 and 1.8 for 
the molecular structures respectively 24J . We therefore 
conclude that at present the diffraction measurements are 
insufficient to distinguish between the proposed dimeric 
and the theoretical molecular PA^/mnm structures. 



Finally, we examine the relative stability of all consid- 
ered crystal structures as a function of both P and T. 
A plot of enthalpies showing that Cmca is never ther- 
modynamically stable at low T is reported in Fig. HJi. 
This finding confirms a previous observation by Iota and 
Yoo that the orthorhombic symmetry obtained dur- 
ing ambient T compression of CO2-I is metastable; the 
reason why the system reverts to it is the low kinetic bar- 
rier associated with the rotation of the molecules from 
their alignment in Pa3 to that in Cmca. In additional 
support of this conclusion, we have performed structural 
optimizations starting from the Pbca structure, and re- 
laxing the atomic coordinates. The Pbca lattice is ob- 
tained from Pa3 by deforming the cubic cell; further ro- 
tations of the molecules in the y-z plane then leads to 
Cmca. In all such calculations above 15 GPa, the CO2 
molecules quickly align in the y-z plane, thus ending up 
in the Cmca, rather than in the energetically favorable 
P^/mnm phase. This happens regardless of whether 
the unit cell is tetragonal or orthorhombic; the deforma- 
tion of P42 / mnm to the orthorhombic Pnnm costs only 
about 1 meV per molecule for 1% distortion of a and b. 

The kinetic barrier between Cmca and P^/mnm can 
be overcome by heating the system above ambient T . In- 
terestingly, we find that upon further heating at constant 
P there is a thermodynamic phase boundary above which 
the orthorhombic phase is now stable. The phonon free 
energies, computed as described in Ref . E3L are plotted in 
Fig- lib- The transition to Cmca is entropy driven and is 
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FIG. 4: (a) Enthalpy versus P relative to the Cmca phase 
at T — 0. (b) Phonon free energy versus T relative to the 
Cmca structure at 16 GPa; the reported differences do not 
vary strongly with P. (c) Relative stability of selected CO2 
structures from Gibbs free energy comparison; here P includes 
phonon contributions. The dotted lines indicate the experi- 
mental phase diagram constraints from Ref. (see Fig. . 
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mainly due to a soft acoustic phonon mode. Its physical 
origin is in the specific geometry of Cmca, allowing for a 
relatively unobstructed shear-like motion parallel to the 
y axis; we have confirmed this by sampling the potential 
surface with frozen-phonon calculations. 

The computed stability regions for Pa3, P^/mnm, 
and Cmca are shown in Fig. Since we are comparing 
the relative stability of a limited number of structures, 
we cannot firmly conclude that Fig.^J: represents the ac- 
tual phase diagram of CO2. However, the almost perfect 
matching of the stability region of Cmca with that of 
CO2-IV makes the identification between the two phases 
rather compelling. We note that a phase IV of CO2 was 
first suggested by Olijnyk and Jephcoat 0] based on the 
observed multiplicity of external Raman-active modes 
between 10 and 20 GPa, at ambient T. Their data is 
consistent with the Raman spectrum of Pa3, Cmca, and 
Pbca. As we discussed above, Pbca is intermediate be- 
tween Pa3 and Cmca, and may be present close to the 
boundary of stability between these two phases, which 
extends from 5 to 16 GPa in our calculations. In addi- 
tion, Cmca has a soft (below 100 cm^ 1 at ~ 20 GPa) 
zone-boundary libron, Bi g , tilting the molecules away 
from the y-z plane, which crosses two of the acoustic 
modes. Thus, when excited at high T, Cmca could actu- 
ally turn into a peculiar distorted phase. R ef . Ilfil observed 
a similar distorted phase with a weak (otherwise forbid- 
den) Raman activity of the bending internal mode, which 
they related to the presence of pressure inhomogeneities 
in the sample. We find such an explanation of the Ra- 
man activity observed by Yoo et al. || around 650 cm -1 
more plausible than the existence of bent molecules. 

In conclusion, we have presented a new interpretation 
of experimental data recently obtained for CO2-II, HI, 
and IV. Our ab inito calculations identify CO2-II as a 
molecular structure with the PA^/mnm symmetry. We 
have also elucidated the high-T behavior of molecular 
CO2, and have found the Cmca symmetry to be prefer- 
able due to its relatively large entropy. Based on our find- 
ings, we propose that measurements of the high-T phase 
IV be interpreted along the lines of an orthorhombic, and 
possibly distorted, phase. Finally, we note that the ex- 
ternal Raman modes of several CO2 molecular structures 
are very close in frequency (see Fig.[3J), and therefore may 
not be sufficient to discriminate between different phases. 
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